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Abstract: We study, in a holographic setup, the effect of adding localized fermionic 
impurities to the three-dimensional Chern-Simons-matter theories with unquenched 
fields in the fundamental representation of the gauge group. The impurities are intro- 
duced as probe D6-branes extending along the radial direction and wrapping a five- 
dimensional submanifold inside a squashed CP 3 . We analyze the straight flux tube 
embeddings and study the corresponding fluctuation modes of the D6-branes. The 
conformal dimensions of the operators dual to such fluctuations depend non-trivially 
on the ratio of the flavor number to the Chern-Simon level of the unquenched ABJM. 
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1 Introduction 

One of the developments emerging from string theory explorations is the idea of a 
gauge/gravity correspondence [1]. The remarkable feature of the correspondence is 
the relation of the strongly coupled regime of the gauge theory to the weakly coupled 
regime of the string theory and vice-versa. Consequently, it has become a powerful 
tool in studying strongly interacting systems, allowing for novel computations that go 
beyond the standard perturbative techniques of quantum field theories. 
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In this theme, the recent developments on the AdS^/CFT^ constitute a rich frame- 
work in which fundamental questions about the correspondence can be posted. In 
particular, the ABJM theory [2] is a U(N) x U(N) Chern-Simons gauge theory with 
levels (k, —k) and bifundamental matter fields. In the large TV limit the theory ac- 
quires a supergravity description in terms of the AdS^ x § 7 /^fc geometry. When the 
Chern-Simons level is large the size of the fiber is small and the system acquires a 
ten-dimensional description in terms of AdS^ x CP 3 with fluxes, that preserves 24 su- 
persymmetries. The addition of flavor to the ABJM theory (fields transforming in the 
fundamental representations (N,l) and (1,N) of the U(N) x U(N) gauge group) is 
realized through D6-branes filling the AdS^ space and wrapping a submanifold inside 
the CP 3 , while preserving a fraction of the initial supersymmetry [3]. The addition of a 
large number of flavor branes, continuously smeared in the transverse space, produces a 
backreaction of the original geometry and induces a deformation. Utilizing techniques 
developed in [4, 5] and reviewed in [6], this unquenched solution was computed in [7] 
and depends non-trivially on the number of flavors. 

The addition of an extra set of branes interacting with the colored ones creates a 
defect in the gauge theory. The characteristic example of this class is J\f = 4 Super- 
Yang-Mills (SYM) with fermionic impurities. The branes realization of this system is 
obtained by adding D5-branes to the AdS§ x S 5 background, dual to M = 4 SYM. The 
D5-branes are extended in the holographic direction and wrap an S 4 inside the S 5 in such 
a way that their worldvolume is an AdS% x S 4 submanifold of AdS& x § 5 [8] (for other 
embeddings with lower dimensional spheres see [9]). More recently, these configurations 
have been used in [10, 11] to construct holographic dimer models, by considering D- 
branes which connect two spatially separated impurities on the boundary of AdS. The 
holographic setup of D3- and D5-branes also realizes the maximally supersymmetric 
Kondo model [12-14] (see also [15, 16]). 

In this paper we consider D6-branes extending along an AdS2 C AdS± and a wrapping 
of a five- dimensional submanifold inside the squashed CP 3 , in order to construct the 
holographic dual of a Chern-Simons-matter theory with flavor and fermionic impurities. 
We will built our solution on the unquenched background solution [7]. 



- 2 - 



An overview of the paper is as follows: In section 2 we give a short, self contained, 
review of the gravity dual of a three-dimensional Chern-Simons-matter theory with 
unquenched fields in the fundamental representation of the gauge group. In section 3 
we introduce impurities as probe D6-branes wrapping a five- dimensional submanifold 
inside the squashed CP 3 and analyze straight flux tube embeddings. In section 4 we 
analyze in detail the fluctuation modes of the D6-branes around the straight flux tube 
configurations and compute the conformal dimensions of the operators dual to such 
fluctuations. Due to the presence on the unquenched ABJM there is an explicit depen- 
dence on the number of flavors. In section 5 we conclude and discuss lines of possible 
future related research. In the appendix A we analytically derive the Lagrangian for 
the fluctuations of the probe brane and in appendix B the explicit solution of the ra- 
dial type differential eq. arising in these fluctuations. In appendix C we compute the 
spectrum of the Laplacian corresponding to the angular part of the operator entering 
into the fluctuation analysis. 

2 Review of the ABJM with unquenched massless flavor 

In this section, following [7], we will review the solutions of type-IIA supergravity 
which are dual to three-dimensional Chern-Simons-matter theories with unquenched 
fields in the fundamental representation of the gauge group. The addition of flavor is 
performed by means of D6-branes that are extended along the Minkowski gauge theory 
directions and are delocalized in the internal space in such a way that the system is 
J\f — 1 supersymmetric and the flavor group is Abelian. The geometry has the form of 
a product space AdS± x M.§, where M.q is a six-dimensional compact manifold whose 
metric is obtained by squashing the Pubini— Study metric of CP 3 [17] with suitable 
constant factors depending on the number of flavors. 

The metric of the flavored ABJM background (in the string frame) is given by 



ds 2 = L 2 ds 2 AdSi + ds 2 & 



(2.1) 



with the standard parametrization for the AdS± metric 



ds 2 A 



2 



-dt 2 + dx 2 + dy 2 ) + 



dr 2 



(2.2) 



AdSi 



= r 



2 ' 



r 
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while the six-dimensional metric is written in terms of the SU (2) instanton of S 4 

L 2 



ds 2 & 



qds 2 & + (dx i + e ijk A jX k ) , (2.3) 



where b and q are constant squashing factors. The metric for the unit S 4 is denoted by 
G?Sg 4 and x % (i = 1, 2, 3) are Cartesian coordinates parametrizing the unit S 2 whereas A 1 , 
i — 1,2,3 are the components of the non-Abelian one-form connection corresponding 
to the SU(2) instanton. The solution depends on two integers N and k which, on the 
gauge theory side, represent the rank of the gauge group and the Chern-Simons level, 
respectively. In string units, the AdS± radius L can be written as 

Ll = 2 ^N (2-,)f (24) 
k q(q + r/q — rj) 

Introducing a set of SU(2) left-invariant one-forms, which satisfy the usual relation 
duji = \ €ijk ojj A ojk, together with the coordinate a (0 ^ a ^ n) one parametrizes the 
unit S 4 as 

ds 2 S 4 = da 2 + Sm a (uf + uj\ + uj\ ) . (2.5) 
The parameterization for the a/s is 

L>j\ = cos ip 76*2 + sin ip sin 82 d<p2 

002 = sin-07# 2 — cos ijj sin # 2 702 (2-6) 
uj 3 = dip + cos 9 2 d(f>2 , 

while the 577(2) instanton one-forms Ai are given by 

Ai = - sin 2 -^. (2.7) 

Parametrizing the x % coordinates of the S 2 by means of the angles 81 and 0i, the 
following relation holds 1 

(dx* + e ijk AjA k ) 2 = E\ + E 2 , (2.8) 
where E\ and 772 are the following one-forms 

Ei = ddi + sin 2 ^ (wi sin 0i — W2cos0i) 

77 2 = sin#i (d(j)i — u^sin 2 — J + sin 2 ^ cos 9i (cji cos 0i + u; 2 sin0i) . (2.9) 



1 Explicitly we have x 1 = sin^i cos^i, x 2 = sin6*i sin^i, a; 3 = cos#i with ^ 6\ < ir and 
^ <j>i < 2w. For completeness we also note that ^ 82 < n, ^ <p2 < 2ir and ^ ip < Att. 



-A- 



Putting all these ingredients together we rewrite the six-dimensional metric (2.3) as 



ds 2 - — 



qdsl 4 + E\ + E\ 



(2-10) 



In order to write the expression for the F 2 we will introduce a new set of one-forms 



S\ = sin 0i oji — cos <f)\ u 2 , 

S 2 = sin 61 U3 — cos 6\ (cos (f>\ u)\ + sin <f>i u 2 ) , 

5*3 = — cos^u^ — sin 6*i (cos 0i Wi + sin^u^) 



(2.11) 



Then, the ansatz for the F 2 is the following 

k 



F 2 = - 
2 L 



E\l\E 2 — 77 (S a A <S 3 + Si A S 2 ) 



(2.12) 



where the one-forms S a and Si are 

S a = da , <Sj 



sin a 



1,2,3 



(2.13) 



and rj is a squashing parameter, directly related to the number Nf of flavors as 

1 , 3iV / 



77 



4fc 



1 ^ 77 < 00 . 



(2.14) 



The internal squashing q is related to 77 through the relation 

3(1 + n) - J9r] 2 - 277 + 9 
9 = 2 ' 

while 6 can be written in terms of g and 77 as 

, g(7? + q) 



(2.15) 



(2.16) 



2(q + rjq — r)) 

The dilaton and the F 4 have the following expressions 

^_b V + qk 3k (77 + g)6 2 

e -4^, 4 ~ T "2^" ^ ^ 4 • (2 - 1?) 

The range of q G [1,5/3], the value q = 1 (no flavors) corresponds to the A/" = 6 ABJM 

background. We also note that had we taken the positive square root in (2.15) it would 

have corresponded to a different branch for which q ^ 5. For the value g = 5 the 

background has reduced supersymmetry, whose metric is the sum of the AdS± metric 

and the squashed CP 3 . The latter metric is not the one corresponding to an Einstein 

space [18, 19]. This would have corresponded to the value q = 2, which nevertheless is 

not allowed in either branch. 
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3 The Hamiltonian density of D6-brane probes 



In this section we will consider a D6-brane probe extending along the radial direction 
r and wrapping a five-dimensional submanifold inside the squashed CP 3 at constant 
values of the spatial Minkowski directions x and y. The background coordinates X M 
and the worldvolume coordinates are 

Background: X M = (t, x, y, r, a, 6> 2 , 2 , ip, #i, 4>i) , 

(3.1) 

Brane: C M = (*, r, 7*) = (*, r , #2, 02, ^, #i, <M • 

We consider embeddings in which the angle a depends only on the radial direction, i.e. 

a = a{r) . (3.2) 

These embeddings correspond to configurations in which the flux tube starts from the 
boundary of the AdS^ and reaches the origin of the holographic coordinate. 2 We will 
also turn on an electric worldvolume gauge field component Fo r , whose source is the 
RR potential C5 through the Wess-Zumino (WZ) term of the D6-brane action. 
The Dirac-Born-Infeld (DBI) part of the D6-brane action is given by 



S 



DBI 



-T 6 J d 7 ( e-V-detQ7 + F) , (3.3) 



where g is the induced metric on the worldvolume of the D6-brane and T 6 is the brane 
tension. After integrating over all the angles of the internal space we arrive at the 
following expression for the DBI contribution to the action 



Sdbi = J dtdr£ DB i, (3.4) 

with 

NL 2 b . 3 / q 



Cdbi = - ^f- sin 3 a J 1 + £ r V 2 - L~ a Fq , (3.5) 
8ir yfq V b z 

where a' denotes da /dr. The WZ part of the action is given by 

£wz = T 6 [c 5 AF= fdtdrC wz . (3.6) 



That excludes hanging flux tubes, namely configurations that reach a minimal value of r and return 
to the boundary. In such cases a non-constant Cartesian coordinate for the embedding is needed. 
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In type-IIA supergravity the RR six-form is defined as the Hodge dual of the RR four- 
form, namely *F 4 = — F 6 . Through that relation we obtain the expression for the 
five-form potential C5 as 

C 5 = - — N C{a) sin Q x sin 2 d6 x A d9 2 A dfa A <# 2 A # , (3.7) 
8 

where 

C(a) = cosa(sin 2 a + 2) - 2. (3.8) 
After integrating over the above five angles 7* we obtain 

N 

Cwz = - — C(a)F 0r . (3.9) 

87T 

The total Lagrangian density is given by the sum of the (3.5) and (3.9) 



C = - — — 

8vr 



L 2 sin 3 a A /l + Jrr 2 a' 2 - L~ 4 F 2 r + ^C(a)F 0r 



(3.10) 



The equation of motion for the gauge field derived from (3.10) implies that 

dC 

■ constant . (3.11) 



dF 0r 

Using the quantization condition derived in [8] one relates this constant to the number 
n of strings (quarks) of the flux tube as 

H = " T " (3 ' 12) 

where Tf is the tension of the fundamental string and n G Z is the fundamental string 
charge carried by the D6-brane. Exploiting the above quantization condition, we obtain 

sin 3 a a/ sin 6 a + C n (a) 2 



1 - L~*Fl + ±ria' 2 V 1 + T? ra 



(3.13) 



where we have defined 

C!Jn,\ = - 

b \ N 
Then from (3.13) we obtain for the field strength 

r 2 a" 2 



C n (a) = ^ (c(a) + ^ I . 



F 0r = L ] + * ' 



y/ sin 6 a + C n 



a) 2 
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In order to eliminate the electric field from the equations of motion we compute the 
Hamiltonian of the system by performing a Legendre transformation in (3.10) 

dC 

U = F 0r —-£. (3.16) 
oF 0r 

Using the above results together with (3.13) we end up with the following formula for 
the Hamiltonian density 



NL 2 b / q 



H = ^ — a/I + i|r 2 a' 2 J sin 6 a + C n {af . (3.17) 

It remains to determine a(r) by integrating the corresponding Euler-Lagrange equa- 
tions. In the next subsection we will constrain our analysis to configurations with 
constant a. Embeddings with a depending on the holographic coordinate are related 
to the baryon vertex of the ABJM theory and will be analyzed in future work (for the 
similar analysis in the AdS§ x S 5 case see [20, 21]). 

3.1 Flux tube configurations 

In this subsection we will calculate the energy density of the configurations with con- 
stant a. Such configurations must satisfy the condition 



m 

da 



= 0. (3.18) 

a'=0 



Since 



with 



&H 

da 



= 3 ElL — sin 3 aA n (a) 



fq 

A n (a) = sin 2 a cos a — ^—C n (a), (3.20) 

the non-trivial configurations with constant a are the solutions of the following algebraic 
equation 

An(«) = 0, (3.21) 
which can be written as a cubic equation in cos a as 

q \ 3 ( 1 3q \ 2q ( 1 In 



j - ( 1 - ^ \ cosa - ^ M - _ ] = 0. (3.22) 
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Due to Bolzano's theorem 3 the function A n (a) has at least one root in the interval 
a G [0, 7r] for every n in the range ^ n ^ N, while the monotonicity of the function 
in this interval tells us that the root is unique. After using (3.20) to express C n (a n ) in 
terms of a n 

C n {a n ) = — sin 2 a n cos a n , (3.24) 
as well as (3.17), we obtain the energy density of the configurations with constant a 



NL 2 b 2 2 / ~ q o , 

E n = — sin a nX cos 2 a n + — sin a n , (3.25) 

Sir q y " 

where a n is a solution of (3. 22). 4 The constant electric field Fq t corresponding to such 

configurations is computed from (3.15). One finds that 

- L 2 cosa n 

for = - (3.26) 
J cos 2 a n + jjz sin 2 a n 

Notice that from (3.22) we have 

A n (a n ) = - Ajv-n(7r - a n ) a N _ n = it - a n , (3.27) 

which combined with (3.25), is telling us that E n is invariant under the change n — y 
N — n, as it should if an object is to transform in the fully anti-symmetric representation 
of the gauge group, with n being the number of boxes in the corresponding Young 
tableaux. The induced metric on the D6-brane worldvolume is 



- r2dt2 + 7^ + ds k 



(3.28) 



which is of the form AdS2 x M 5 , with the line element of M 5 having the following 
expression 

ds 2 M5 = gijdfd^ = ^ sm 2 a n (u\ + co 2 2 + u 2 ) + i (E 2 + E 2 ) . (3.29) 

Under the change of a n described in (3.27) the line element of (3.29) remains invariant 
as well. 



3 Note that q/b 2 G [l, j|] , monotonously as 77 € [1, 00) and that 

A " (0) = < ' AnW = ^ i 1 ~ v) > ■ (3 ' 23) 

4 For n«J\f the energy turns out to be simply the sum of the energies of the individual fundamental 
strings, i.e. E n ~ r^L 2 , where the extra factor L 2 is related to the overall appearance of the same 
factor in (2.1) due to a rescaling of the world-volume coordinates. It is worth noting that in this dilute 
type approximation there is no dependence of the energy on the flavor number. 
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4 Fluctuations of the impurities 

Moving one step forward, we will study in this section fluctuations around the static 
configurations we have computed. For this reason we consider the following 

a = a n + £(C) , F = f + /(C) , x = x + X (C) , (4-1) 

where a n is a solution of the condition A n (a) = 0, x is the constant Cartesian coordinate 
of the unperturbed D6-brane and / is the background gauge field strength with non- 
vanishing component given by (3.26). The fluctuations around these constant values, 
namely £, x an d /, depend, as indicated, on the D-brane coordinates C M i n (3.1). 
The total perturbed Lagrangian density is the sum of the DBI and the WZ parts, and 
a detailed derivation is presented in the appendix A. Indeed, if in (A. 23) we neglect 
constant and total derivative terms we find the following Lagrangian density for the 
quadratic fluctuations 

c = -n L2 p 1 / 2 { \l 2 r 2 d, x + l 2 g^t 

+ 1 _g» P g»« f ^ fpa +Ve-WUor) , (4.2) 

where we have for notational convenience set 

n sin 6 a n 3 



sin 6 a n + C n (a n ) 2 ' 2 sin 2 a n 



W = r9 m/ 9 -3 — C«(«n) cota n + — sin a n ) (4.3) 
L 2 P 1 ' 2 sm d a n V b ) 

and the seven-dimensional metric Q is defined in (A. 7). 
4.1 Fluctuation of the Cartesian coordinate 

Here we study the fluctuations of the Cartesian coordinate which do not couple to those 
of the gauge field and of the embedding angular coordinate. The equation of motion 
for x can be easily derived from (4.2) and it is 

d r {r%x) ~ dlx + r 2 P V\ h x = • (4.4) 

To obtain the above equation of motion we explicitly used the components of Q^ u , while 
V 2 ^ is the Laplacian operator of the five-dimensional manifold in (3.29) (see appendix 
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C). The actual expression for this operator is quite complicated. Remarkably, we were 
able to find explicit solutions of the form either \ — x(tj r ; 82, 02, ^) or X — X^i r , 4>i) 
(see appendices C.l and C.2, respectively). Without loss of generality in the following 
analysis of the conformal dimensions we will focus on the first class of solutions. Using 
the separation of variables 

X = e* Et R(r)n(0 2 ,4> 2 ,i>), (4.5) 
and equation (C.5) from appendix C.l 



b 2 b 2 1(1 + 2) 
V^ g Q = — — — V 5 3^ = t^- 1 fi, I = 0, 1, 2, . . . (4.6) 



the equation of motion for the radial function R(r) becomes 

d r (r*d r R) + ( E 2 - b l l l±^lp r A R = o. (4.7) 
V <? sin z a n J 

This equation can be solved exactly, but since we are interested in the asymptotic 
behavior of the solution we put all the details on the analytic derivation in appendix 
B. Assuming that R(r) ~ r A at large r, we arrive to the following quadratic equation 

A(A + 3) = ^^P, (4.8) 
q sin a n 

with solutions that we will denote as Ai and A2. We would like to associate them with 
the dimensions A of operators on the defect theory. Since the fluctuations are not 
canonically normalized since there is a factor r 2 in front of the kinetic term for the field 
X in (4.2). Hence, we cannot simply use the usual relation [23] 



with d — 1. Instead one may employ an approach that gives the result immediately 
[24] According to this, if a scalar field in AdSz at large r behaves as 



X 



d ir ~ 2Xl + d 2 r~ 2X \ A 2 >Ai, (4.10) 



then the dimension of the operator dual to the normalizable mode is 

A = ^ + A 2 -A 1 . (4.11) 
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In our case the conformal dimension becomes 




* + p _l , 9 | 1(1 + 2) 

q sin 2 a n 2 V 4 q/b 2 sin 2 a n + cos 2 a r , 

In general, the conformal dimension is not a rational number and A = 2 for I = 0. 

Moreover the dimension depends on the filling fraction u, through the dependence of the 

angle a n , and is invariant under the change of (3.27). Also the range of the parameter 

q/b 2 in footnote 2, the prefactor multiplying 1(1 + 2) is of order one. 5 The dependence 

of the conformal dimension on the number of flavors becomes is in general complicated 

and it admits a power series expansion around the unquenched result. In particular, in 

the half filling fraction case, v — | the cubic equation (3.22) has the following solution 

in the interval a n G [0, 7r] 



7T 



cosa„ = a n = - . (4-13) 



Expanding (4.12) around the unquenched result we have 



1 7 9 , ui , on 9 1(1 + 2) f N f\ „ ( N f 



A = — h \ - + 1(1 + 2) , -r- +0[-±) . (4.14) 

2 V4 512^9 +i(/ + 2) 

4.2 Coupled modes 

In this subsection we will focus our attention on the fluctuations of the gauge and scalar 
fields, which through their equations of motion appear to be coupled. The equation of 
motion for the gauge field is given by 

d p ( G^G^Uu) + W ( d r t 5° - d Z 6?) = 0, (4.15) 

while that for the scalar is 

^8,G/S<T«) - ^Vi + £±Wf» = 0. (4.16) 

We consider the following ansatz for the fluctuations of the gauge fields and the scalar 
(the only non-vanishing components of the gauge field are A r and Aj) 

A r = e lEt Q(6 2 , 2 , iP) <f>(r) , ii = e lEt dtf^ 02, ^) 4>(r) 

^ = e^fi(^ 2 ,0 2 ,^)z(r). (4.17) 



5 Another way to arrive at the same result is to redefine the fluctuations by absorbing the factor r 
in the kinetic for x i n (4-2) into a new scalar field 4> = r\. Then <f> becomes canonically normalized 
but after some algebraic manipulations one sees that m 2 is shifted by a factor of 2. Then after using 
(4.9), with d = l, one derives (4.12). 
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Note that the fact that the vector index of Aj is due to the derivative on Q is, as it 
turns out, for consistency of the ansatz upon substitution into the above equations of 
motion. Then, from (4.15) for a = we obtain 

2 
P 



b 2 l(l + 2) ~ 



— L 4 PWr 2 z, 
E 



q sin a n 

while if we set a = r in the same equation we have that 



(4.18) 



E 2 M+ b 2 l(l + 2) r 2 



L*P 2 ^ ' q sm 2 a n L 4 P 
Also from (4.15) for a = i we have 

o d 



iEW z = 



E 1 



dr 



r 2 [(f)' 



(4.19) 



(4.20) 



an equation that can be easily derived from (4.18) and(4.19) . Using (4.18) in order to 
eliminate <fi from (4.19), we have 

,dz\ E 2 



dr 



dr 



r z q sin a n E 



d 
dr 



dr 



+ 



while the equation of motion for the scalar, using (4.17), becomes 



d 



dz\ E 2 



dr \ dr J 



+ 



<1 



2V 



1(1 + 2] 
sin 2 a„ 



b 2 

P z + i—E PW(j) 



0, (4.2i; 



(4.22) 



At this point, we define the differential operator O, which acts on functions as follows 



d 



df 



dr \ dr 

and make the following field redefinitions 



E 2 



L 4 

z = i—P 2 W z 
E 



r ] = ( f ) -i P 2 Wz 
E 



— z . 



Then (4.22) and (4.25) can be written in a more compact form as 



CI 2 x2 - M 



(4.23) 



(4.24) 



(4.25) 



where the entries of the matrix Ai are 

q 
b 2 



2V + L i W 2 P 2 + /(/ + 2) 



sin 2 a r 



M 



12 



q 



-L 4 W 2 P 3 , 



M 



21 



M 



22 



b 2 1(1 + 2) 
q sin 2 a n 



(4.26) 



13 



while its eigenvalues are 
P b 2 



2V + L 4 P 2 W 2 + 2 ± J4 l -V±^-L±P 2 W 2 + (2V + L 4 P¥ 2 ) 



sin a n y sm a n 

(4.27) 

Now, if are the eigenvectors of yVf with eigenvalues A±, then (4.25) takes the form 



dr V dr / V r 



r 2 ^) + (^-A ± )^ = 0, (4.28) 



and in order to study the behavior of ip± at large r we assume that ip± ~ r s . Then 
from the above differential equation we find that s should obey the following quadratic 
equation 

s(s + 1) - A± = 0, (4.29) 

with solutions 

-1 ± Jl + 4A± 
s± = \ , s_ < s + . (4.30) 

Noting that in the kinetic term in (4.2) for these type of fluctuations there is no overall 
extra overall factor of r, we may safely use for the conformal dimension the expression 
[23] with d — 1 and m 2 = A+, i.e. 

A = X - (l + Vl + 4A+) . (4.31) 

In general, it is not a rational number and depends on the filling fraction v. Notice, that 
unlike (4.12), for I = the conformal dimension A + does depend on the flavor number. 
As also argued in the appendix C, for angular dependence of the form f2(#i, 0i) we may 
use the previous results by simply performing the replacement (C.8). In the half filling 
fraction case v = | and expanding (4.31) around the unquenched result we obtain that 

9 (/ + 3)(( + 2)(2(-l) fNA 2 , „(Nf- 



A = 3 + < - 512 (/+W+20 [t) + °{t) < 432) 
5 Conclusions and future directions 

In this paper we studied localized fermionic impurities to the unquenched ABJM 
Chern-Simons-matter theory, which are realized through the addition of fields trans- 
forming in the fundamental representations (N, 1) and (1, N) of the U (N) x U (N) gauge 
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group. In the holographic approach the impurities are added by introducing probe D6- 
branes, extending along the holographic coordinate and wrapping a five-dimensional 
submanifold inside a squashed CP 3 at constant values of the Minkowski directions. The 
background RR field induces an electric gauge-field on the world-volume of the probe 
branes, giving rise to a bundle of strings that form a flux tube which prevents the 
collapse of the wrapped brane. 

We concentrated on the simplest of the configurations in which the flux tube starts 
from the boundary of the AdS± and reach the origin of the holographic coordinate. 
More general solutions including the baryon vertex of the unquenched AB JM theory and 
hanging flux tubes, namely configurations that reach a minimal value in the holographic 
coordinate and return to the boundary, are left as open problems for future work. 

The natural step forward was to investigate of the stability for the probe D6-branes, 
that introduce the holographic impurities. We presented an analytic study for the 
fluctuations of those probes in the unquenched ABJM. The fluctuations are separated 
in two categories. The first contains just the decoupled fluctuations of the Cartesian 
coordinates while the second the coupled fluctuations of the angular embedding function 
and the world-volume gauge field. The coupled modes were whown to decouple by 
appropriate field redefinitions. In this way we were able to determine the spectrum of 
conformal dimensions of the dual operators in the defect theory. The novel feature of 
our analysis is that using the unquenched ABJM background we obtained expressions of 
the conformal dimension that explicitly depend on the number of fundamental flavors, 
thus generalizing the previously obtained results for the ABJM background [11]. 

There are many interesting questions that follow from the analysis we presented 
in the recent paper, and we would like to pursue some of them in the near future. In 
the quenched ABJM background there is robust proposition in [22] that the Wilson 
lines in the antisymmetric representations of the gauge group of the ABJM are holo- 
graphically realized by D6-branes extending along an AdS2 C AdS^ and wrapping a 
five- dimensional submanifold of CP 3 . These are the natural brane configurations to 
consider in order to construct the holographic dual of the Chern-Simons-matter theory 
with fermionic impurities. Unlike to the quenched case, there is no such proof for the 
unquenched ABJM, though we note the invariance under n — > N — n in (3.27). It 
would be very interesting to pursue this issue further. 
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The addition of fermionic impurities in the unquenched AB JM background at finite 
temperature [25] will create a much richer structure. The analysis of the thermody- 
namic properties, of both straight and hanging flux tubes, is expected to unveil a 
competition between the two configurations. This in turn will lead to the existence of a 
dimerization transition similar to the one presented in [10, 11], but now in a background 
that will include the non-trivial presence of fundamental flavors. 
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A Fluctuation analysis 

In this appendix we will analyze the small perturbations around the flux tube configu- 
rations, derived in section 3.1. We will analytically obtain the second order lagrangian 
for those fluctuations, which is the starting point of section 4. 

We perturb a D6-brane probe as in (4.1) and expand the DBI+WZ action to second 
order in the perturbations £, / and x- Starting with the components of the perturbed 
induced metric we write 



where g is the zeroth order induced metric and the perturbation g has the form 



9 = 9 + 9, 



(A.l) 



g, u = L 2 \r 2 d, X duX + jz d £ d »Z + ^ + ? 9$] +••• , 



(A.2) 
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where the g$ and g$ are given by (their indices take values only in the angular part) 



i(2) 



9if dy = -j- sin 2a n ( u{ + u 2 2 + w 2 ) + 



Q_ 

Ab 2 



b 2 



_ <9Fi _ dE 2 
E\— — + E 2 



da 



da 



9$ *f d^ 



Ab 2 
1 

b 2 



cos 2a n ( cof + + ) 



(A.3) 



9a 



9a 



<9 2 E 1 



+ F 



<9 2 F 2 



9a 2 da 2 

The determinant of the DBI part can be written as 

det(s + F) = det(s + /) det(l + X) , X = (g + f )- l (g + /). (A.4) 

Hence the important step is to compute the components of matrix X in the expansion 

V / det(l + X) = 1 + - TrX - - TrX 2 + - (TrX) 2 + 0(X 3 ) . (A.5) 

2 4 8 

The matrix ( g + / )~ x can be written in a block diagonal form 



(9 + f)- 1 



Q- 1 + J\ 



Or 







1'J 



(A.6) 



where Q 1 and are its symmetric and antisymmetric parts, respectively. The non- 
zero elements of those matrices are 



sin 6 q w + C n (a r 
r 2 L 2 sin 6 a„ 



and 



J 



Or 



J 



rO 



sin 6 a n + C n (a n ) 2 2 

T" 

L z sin a„ 



C n (a n ) ^sm & a n + C n (a n )' 



(/'•' = /. 2 //'-' 



L 2 sin 6 a n 

The matrix elements that contribute to the TrX are 
sin 6 a n + C n (a n ) 2 



(A.7) 
(A.8) 



A o 



r z sin a r 



r 2 (d oX y 



IP 



my 



+ 



C n (a n ) Vsin 6 a„ + C n (a n ) 2 



sin 6 a n 



r 2 d x d r x + p <9 { d r i - L 2 f 0r } 



XI 



sin 6 a n + C n (a r 



sin a r 



r 2 { r 2 («9 rX ) 2 + ^ (9,0 



(A.9) 



Cn("n) V sin «n + C 'n( 



sin 6 a 



I r 2 <9 x d r x + p <9 £ <9 r £ + L 2 f 0r J 



Xj = r 2 ^4 X ^X + J^Qtffirf + £ (MW); + e 2 (M( 2 ))J + L- 2 ^ fe /^ 
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where we have defined the matrices 



(M«);. = y*$g, 2 = 1,2. (A.10) 

The metric cjij corresponds to the five-dimensional space whose line element is given 
by (3.29). In order to calculate the trace of X 2 , we need to compute the non-diagonal 
elements of X up to first order in the fluctuations. The matrix elements that contribute 
to TrX 2 are 

_ _ sin 6 q w + C n (a n ) 2 
r A L l sin a n 

vr sin 6 a n + CJ 
L z sin a„ 



sin 6 a n + C n (a n ) 2 C n (a n ) ^ sixf a n + C n (a n ) 2 

A i o T1 ■ 6 JOi H r2 • 6 l A - ii J 

1/ sm a n L z sm a n 



XI = L~ 2 ~gVf j0 

_ sin 6 q n + C n (a n ) 2 2 C ra (a n .) y / sin (i a w + C n (a n ) 2 

■X-i — T o ■ 6 r Jri r2 • 6 

L z sm a n L z sm a n 

XI = L~ 2 ~g^f jr 

Putting everything together we find that the TrX is given by 



m v n C n (a n ) ^/sin 6 a n + C n (a n ) 2 2 1 + (3g - l)cos2a n 2 

TrA = —2 L / 0r + 6cota„4 H ^ £ 

sin a n g sin a n 

+ L 2 r 2 ^^ X ^x + J^S^e&e, (A. 12) 

where we have used the following equations 

TrM« = 6 cot cu , TrM (2 ^ = \±^Z^^ . 



(A.13) 
q sm a n 



The Tr(A 2 ) is given by 



rp (v2 s _ or _ 4 sin 6 a n + 2C n (a n ) 2 sin 6 a n + C n (a n ) 2 2 2 1 + 3g + (3g - 1) cos 2o^ 



sm a n sm a n g sm a n 

+ L - 4 i 2 sin 6 a n + C n (a n ) 2 fl _ ^ sin 6 a n + C n (a n ) 2 ^ ^ _ ^ \ ^ M ^ 

sin 6 a„ r 2 sin 6 a„ ™ fcj ' J 
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while for the (TrX) we have 



frn v\2 at -4 Snl a n ^n( ft n) ri I \2 r2 , oc j-2 tl 

(TrX) = 4L ^ C n {a n ) fa + 36 cot a n £ 



sin a 



r> A T -2 C n 

(a n ) Jsm & a n + 

- 24L r-g cota n 4 f 0r . (A. 15) 

sin a n 

Putting everything together we calculate the expression of y/det(l + X) 



n T7Z , y x , C„(a n ) ^/sin 6 a n + C n (a n ) 2 2 1 2 2 

Vdet(l + X) = 1 * — _^ L / 0r + 3 cot a n £ + - L r Q 11 d^vX 



sin a 



+ \± I?Q^d£d v t + ? (2 cot 2 a n - 1)£ 2 + ~ Q^QTMp* (A.16) 

o r -2 C n (a n )^J sm & a n + C n (a n ) 2 

- 3L r-g cot a n 4 /or + . . . . 

sin a n 

Since 

det(g + f) = - . r Tr a ( v ■ det ~9, (A.17) 
sm D a n + C n (a n y 

the DBI part of the Lagrangian density is 

_ 7r 2 N b 6 L 2 sin 3 a n r- n — ; — - 

C DBI = -T e - 2 == - vWdet(l + X), (A.18) 

Q V Sln «» + ^n(«n) 2 

where g is the determinant of (3.29) 

r <1 3/2 

\/g = — — sin 3 a n sin#i sin 6*2- (A. 19) 

8 o 5 

What remains is the computation of the WZ part. Using the conventions of section 3 
we have 

7T 2 X 

C-wz = -T 6 — - sin^ sin 6 2 C{a) F 0r , (A.20) 

o 

where F 0r = f 0r + f 0r and the function C(a) has to be expanded around a n 

9 

C(a) = C(a n ) — 3sin 3 a n , £ — - sin 2 a n cosa„ £ 2 + . . . . (A. 21) 

Putting everything together, the WZ part becomes 

C-wz = ~n^- sin 6 1 sin fl 2 1 L 2 C w (a w ) ^ ^ + ^ (A ^ 

° [ V sm "n + C n (a n ) 2 

3L 2 sin 3 a n C n (a n ) 3 9 2 sin 2 a n cosa n C n (a n ) 

£ - 3 sm a n / 0r £ - -L ===== £ 



V^in 6 a„ + C n (a n ) 2 2 A/sin 6 a„ + C n (a n ) 2 
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Finally, summing the DBI and WZ parts we obtain the result 



7T 



WL 2 sin 3 a n , y/qC n {a n 



C DB i + C W z - -T 6 g / . fi ^ / x V 9 i 1 + ~T — r-g C(a;„) 

<T Vsm 6 a n + C n (a n ) { b sm b a n 



b 



sin av 



+ 



( |(2cot 2 a n - 1) - ^f^ n (a„)) £ 2 + V^/,^ 
\2 2 6 sm a n J 4 



3L — g- ( sm a n + C n (a n ) cot a n ) f / 0r f . 



y/sin 6 a n + C w (a n ) / ^ 
sin 6 a„ \ b 

In the above action the first two lines can be dropped containing either constant or 
linear or total derivative terms. The term linear in £ also vanishes upon using the 
condition (3.21). The remaining terms constitute the action (4.2) used in the main 
text. 

B Solution of the differential equation for R(r) 

The differential equation (4.7) is of the form 

T-( r4 rr) + (E 2 -Ar 2 )R = 0, (B.l) 
ar \ dr J 

where A is a constant. To solve this differential equation we first study the asymptotic 
behavior of R(r) at large r. Setting R(r) ~ r A we end up with the condition 



X(X + 3) = A =► X ± = - 3 -±^ 9 - + A. (B.2) 

Now that we know the asymptotic behavior of R(r), we make the ansatz R(r) = r A /(^) 
and obtain the following differential equation for the function / = f(u) = /(1/r) 

u 2 P, - 2{\ + l)u^- + E 2 u 2 f = 0. (B.3) 
du A du 

Now we observe that if we substitute f(u) = u x+ ^g(u) the above equation becomes 

u 2 p- 2 + u^- + {E 2 u 2 -t 2 )g = 0, (B.4) 
du z du 
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where t 2 = (A + |) 2 = § + A , which is the Bessel differential equation with solution 

9{u) = CrJ^iEu) + C 2 Y^(Eu) . (B.5) 
Finally, the solution for R(r) is 

IHr) = (f ) + fty^j (f )) , (B.6) 

which has the correct asymptotic behavior at large r 

R(r) = Cir-i + v / !+ I + C 2 r~ f ~ v/ ?+ I . (B.7) 
In our case the constant A appearing in (B.l) is 

A = b l'lpl P . (B.8) 
q sin a n 

C The Laplacian on the five dimensional manifold M5 

The action of the Laplacian of the five-dimensional manifold M 5 , with metric g^ (3.29) 
on a scalar depending on all the coordinates of M5 is 

v 2 m 6 / = r^krVla/ + °-f-^d ei (sinft^/) + c -f-\rdU 

4 sin ^ 2 sin ft 2 sin ft 

_ £l£ in^- i ) _ 

sin ft sm ft 

Ci cos(0i - ip) Ci cos(0i - V) „ o f. 
- sin ft sin ft ^ sMft— (C ' 1} 

Cl COs(0i '0) COt ft o , . a r, ,x . 4./J q / /j. no « 

H — a o dl (smftcfy/J + ci cot ftcfy (cos(0i - ip)d dl f) 

sm ft 

j. /i o / /. I\X f\ Ci cos(0i - ^) cot ft . 

-ci cot ft (cos(0i - ip)o 02 f) — d d2 (smftc^/) 

sm c?2 

, cot #i<, / ■ / / (U cot ft sin(0i - 0) 

+ci— (sm(0i - rl>)dfaf) + c x — o^d^J 

sm (72 sm t>2 

+ . /1 Cl . „ dfa ((sin ft sin ft - cos ft cos ft sin(0i - -0)) d^f) 
sm ft sm ft 

/) Cl - n ^ (( sin ^i sin #2 - cos ft cos ft sin(0i - ip)) dfof) 
1 ft sm ft 



+ — 

sml 
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where the last term is the Laplacian of the 3-sphere 

V|s/ = 4 \^-de 2 ( S me 2 de 2 f) + ( -\d^ - cotfcfy) f + d\] 
I sin & 2 \ sm "2 J 

and the constants Ci, c 2 are given by 

b 2 1 

Cl = 7 cos 2^ ' c 2 = 1 + g + (g - 1) cosa n . (C.3) 

Solving the corresponding eigenvalues problem in full generality seems very difficult. 
Instead, we look for configurations of the function /, in which the Laplacian acts on, 
that do not depend on all variables. We found two consistent truncations that simplify 
our eigenvalue problem and we present them in the following subsections. 

C.l Solutions of the form / = f2(# 2 ,02,V ; ) 

In the analysis that we presented in the previous sections of this paper we focused 
on configurations of the form / = f2(#2, 02, V0- For such configurations the above 
expression for the Laplacian becomes 



(C.4) 



CiCOs(0i-^) CiCOs(0i -1p)cot9 2 r ) ( . n a n \ 

sm U\ sm O2 sm v\ 

- ci cot Oidfa (cos (0i - 4>)de 2 ty + c x CQt ^ (sin(0i - ?/>)<9 02 fi) 

sm 02 

+ ~ — — TT^i (( sin ^1 sin ^2 - cos 6 1 ! cos 02 sin(0i - ip)) d^Q) . 
sm 6\ sm 02 

Notice that there are still derivatives to be taken with respect to the angles 61 and 0i 
as well as explicit dependence on these angles. Consistency of our truncation requires 
that all such dependencies drop out completely which indeed turns out to be the case 
as we are left with 

v^ fi = r^kr • ( c - 5 ) 

4 sin ! y 

In this way our eigenvalue problem is reduced to that of V| 3 , a well known operator, 
which has eigenvalues A = —1(1 + 2) where I — 0, 1, 2, . . .. 
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C.2 Solutions of the form / = f2(0i,0i) 

Another truncation that simplifies the eigenvalue problem is an ansatz of the form 
/ = fi(6 l 1 ,0 1 ). Then (Cl) simplifies as 



7 2 n C 1 C 2 1 o / • Q o n x . ClC 2 1 2 Cl sm(0! - ip) 

5 2 sin 0i 2 snr 9 1 Vl sin 2 



ci cos(0i — -0) cot 



sin 6 1 



-e\, (sin 02001 ^) + ci cot 2 5^ (cos(0i - ^deS 1 ) (C6) 



2 



H — : — — 7p<9w, ((sin 0i sin 2 - cos 0i cos 2 sin(0i - ?/>)) ckfi) . 
sin 0i sin 2 

Again all dependence on 0i, <f>\ and ip, finally obtaining that 

= IT ST** (8in + ^S^" = IT v - fi ' < C7 > 

We see that configurations of the form / = /(0i, </>i) lead to the well known eigenvalue 
problem of the Laplace operator on the unit S 2 , with eigenvalues A = —1(1 + 1), where 
1 = 0,1,2,... 

The results we have obtained for the fluctuations using the ansatz / = 0(02, 02, 
can be trivially extended to the case when / = O(0i, 0i). We simply have to compare 
(C.5) and (C.7) and make the replacement 

1(1 + 2) 4 sin 2 ^ (sin 2 ^ + gcos 2 ^) Z(Z + 1) . (C.8) 
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